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We study the double exchange model on two lattice sites with one conduction electron in the 
limit of an infinite Hund's interaction. While this simple problem is exactly solvable, we present 
an approximate solution which is valid in the limit of large core spins. This solution is obtained by 
integrating out charge degrees of freedom. The effective action of two core spins obtained in the 
result of such an integration resembles the action of two fractional spins. We show that the action 
obtained via naive gradient expansion is inconsistent. However, a "non-perturbative" treatment 
leads to an extra term in the effective action which fixes this inconsistency. The obtained "Berry 
phase term" is geometric in nature. It arises from a geometric constraint on a target space imposed 
by an adiabatic approximation. 



I. INTRODUCTION 

A double exchange model was introduced in |Q, |j to 
describe the motion of conduction electrons in the back- 
ground of magnetic ions. In its simplest form the Hamil- 
tonian for the model can be written as: 



H 



h.c.) - y^^JHStSj. 



(1) 



where c; is the electron annihilation operator on the site 



of the lattice, Si 



\(JCi is a spin 1/2 of an electron 



|l2| at site i, Si is a. core spin. The first term of 
describes the hopping of conduction electrons between 
nearest neighbor sites i and j of the lattice. The hop- 
ping amplitude is t. The second term is a ferromagnetic 
Hund's interaction between the core spin and the spin of 
the conduction electron at the same lattice site. To spec- 
ify the model completely we have to fix the total number 
of conduction electrons in the sample. We denote filling 
factor X = Vl- c]ci/N , where A^ is a total number of sites 
of the lattice. In this paper we are interested in the 
limit of a very strong Hund's interaction Ju/i oo and 
of very large core spins 5 S> 1. In (0) we did not include 
any (bare) core spin interactions, e.g., due to superex- 
change. The only interaction between core spins in the 
model occurs through the exchange by conduction 
electrons. 

Because of an infinite Hund's coupling Jh an electron 
spin is always aligned with a core spin on any given lat- 
tice site. Therefore, the effective hopping amplitude of 
an electron between two nearest neighbors is the largest 
when corresponding core spins are parallel. In this case 
an electron can hop between lattice sites without chang- 
ing its spin orientation. One concludes that conduction 
electrons induce a ferromagnetic interaction between core 
spins d, 1. 

In the limit of large core spins 5 3> 1 we consider one 
can use an adiabatic approximation and think of core 
spins as of slow variables so that the conduction elec- 
trons (fast variables) are always in the ground state in 



the background of slowly moving core spins. It is tempt- 
ing to average over the motion of conduction electrons 
(to integrate conduction electrons out) and obtain an ef- 
fective spin model describing the dynamics of core spins 
(and electron spins) induced by charge degrees of free- 
dom p^ . In an almost ferromagnetic state the average 
spin per lattice site is S + corresponding to the sum 
of a core spin and an average electron spin. Therefore, 
heuristically, one arrives to the model of ferromagneti- 
cally interacting fractional spins. However, this immedi- 
ately raises the question of what one means by a frac- 
tional spin. Although the average spin per lattice site 
is fractional, the total spin operator on a given site may 
take only values S (no electron on the site) or S + 1/2 
(there is an electron on the site). Moreover, a fractional 
value of 2S is incompatible with quantum mechanics of 
spins (see section ||). Resolving this apparent contro- 
versy is a main goal of this paper. 



The paper is organized in the following way. The aim 
of the section || is mainly pedagogical. We establish no- 
tations and remind the reader how to write the path in- 
tegral for a single spin. We show how the ambiguity of 
the Berry phase for a single spin requires 25* be an in- 
teger. For completeness we derive the Berry phase for 
a single spin in Appendix ^ using a fermionic determi- 
nant representation for the action. The read er fa miliar 
with the subject can start from the section |ll| where 
we consider the double exchange model (|l|) in the limit 
Jn/t —>■ oo and S ^ I. In section ^ we formulate a 
double exchange model on two core spins and obtain our 
main result - an effective action for core spins induced 
by a single conduction electron. This action (namely, its 
kinetic part) answers the question of how to write Berry 
phase for a system of (two) fractional spins. Appendix ^ 
contains more accurate calculation of the effective action 
which uses the method of gradient expansion. We discuss 
the topology behind Berry phase for fractional spins in 
section ^ and conclude in section 
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II. PATH INTEGRAL FOR SINGLE SPIN 

In this section we remind the reader how to write down 
the path integral for a single quantum spin and discuss 
topological reasons for spin quantization. 

Quantum spin S is defined as a representation of di- 
mension 25' + 1 of an SU{2) Lie algebra with commuta- 
tion relations 



2S\ 



(2) 



where as usual 5^ = 5^ ± iSy. An SU{2) Lie alge- 
bra has irreducible finite dimensional representations of 
dimension 25* + 1 labelled by spin S taking integer or 
half-integer values so that — S{S + 1). We would 
like to construct the classical action for the time depen- 
dent unit vector n{t) which upon quantization by path 
integral reproduces (H). One could start directly from 
and construct the path integral using coherent states 
corresponding to an SU{2) group (see e.g., Ref.Q). In 
Appendix ^ we give an alternative derivation of this path 
integral using the fermionic determinant representation. 
The result is given by partition function 



Z = 



Dn e 



-We 



(3) 



where 



WE[n] = i{2S)2T: [ cf x —e^^'^e^^n'' 8^,71^8^71" . (4) 
Jb 

Let us explain what (^||) mean. First of all we used 
an Euclidean time representation replacing real physical 
time t by an imaginary time r = it. The change to an 
imaginary time changes quantum amplitude e*^ used in 
path integral into a "Boltzmann" factor e"'^^. We im- 
pose periodic boundary conditions (in imaginary time) 
on n(r) so that n{[3) — nifi) and /? is a maximal span 
of an imaginary time. With these boundary conditions 
one can think of time as of a one-dimensional circle with 
identification of point t = [3 with t — Q. Physical results 
will be obtained in the limit /3 — s- 00 after an analytical 
continuation back to the real time t = —it. Imaginary 
time representation is slightly more convenient for calcu- 
lations and for a discussion of topological aspects of spin 
quantization. 

Secondly, the two-dimensional integration in (^) is per- 
formed over some closed domain B of two dimensional 
plane such that the boundary of this domain is identi- 
fied with an imaginary time. The values of n field in 
the area B are defined as an extension from the physi- 
cal trajectory n(r) so that n(r, u) is a smooth function 
of T G [0,/3] and auxiliary variable u E [0, 1] such that 
n{T,u = 0) = (0,0,1) and n(r, u = 1) = ^(t). We 
shall refer to the action (^) as to a Berry phase B for 
a spin (see Appendix ^). Although the action (Q) is 
given in terms of two-dimensional integral over auxiliary 



space, the variation of the integrand is a full derivative. 
Therefore, the variation of Berry phase (^) is given by a 
one-dimensional integral over the "physical" (imaginary) 
time as 



5We = iS I dTt'''"'5n''n'drn' 



(5) 



We see that the variation of (Q) does not depend on the 
choice of an extension of n(r) in the unphysical domain B. 
Using (|^) one can show that the equation 5{We + W^) = 
0, where = dr Shn is a coupling to an external 
magnetic field gives a classical equation of motion ic^T-n = 
[h X n]. The latter is the equation of precession of spin 
in magnetic field h. 

Partition function (||) with the action (Q) solve the 
problem of path integral representation for the algebra 
(^. Namely, the commutation relations (|^) can be ob- 
tained as a result of quantization of "classical action" 



A. Berry phase and quantization of spin 

We have already checked that an infinitesimal variation 
of Berry phase 5We [n] does not depend on the extension 
of ^(t) into the auxiliary domain B. How about VF£;[ri] 
itself? One can easily show (see Appendix ^) that the 
action (Q) is equal to i{2S)^/2 where J7 is a solid angle 
swept by a unit vector n during its motion (see Fig. |^). 
This angle is not uniquely defined: one can always add 
or subtract any angle which is a multiple of a full solid 
angle Att. 

Therefore, W£;[n] is not unambiguously defined and 
is a multi- valued functional of ^(t). Multi- valued func- 
tionals of this type were studied in both mathematics 
and quantum field theory (see e.g., [^) and are often 
referred to as WZNW term — Wess-Zumino-Novikov- 
Witten term. Although WE{n) is a multiple- valued 
functional, in partition function (^) one needs e~^^["l 
which is a single-valued functional of n(T) under the 
condition that the coefficient 2S is an integer. Then 
^~5We ^ g-i(2S)27rfc ^ ^ ^^^^ Boltzmauu weights 
Q~WE{n) defined unambiguously by physical configu- 
rations n(T). One can think of this constraint on the 
value of 25* as of the reason for spin to be an integer or 
a half-integer. Having non-integer 2S in (^) amounts to 
having inconsistent theory. 



B. CP^ representation of Berry phase and topology 

There is a one more representation of Berry phase 
which is particularly convenient for the purposes of this 
paper. It is a so-called CP^ representation. We write n = 
zWz, where z = {zi^z^f is a complex two-component 



vector satisfying the constraint z^z 
The latter reproduces = 1. 
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n(T,u=0)=(0,0,l) 




or even more formally that the second cohomology group 
of 5^ is H2{S^) = Z. This allows one to write down 
the topological term where an integration of the men- 
tioned 2-form is performed over a domain B with physical 
time as the boundary of this disk. The special proper- 
ties mentioned above guarantee that (Q) unambiguously 
defines the path integral (0). 




FIG. 1: The unit vector n(T) draws a closed line on the sur- 
face of a sphere with unit radius during its motion in imag- 
inary time. Berry phase is proportional to the solid angle 
(shaded region) swept by the vector n(r). One can calcu- 
late this solid angle by extending n into the two-dimensional 
domain B as n{u,T) and calculating (^). 



In terms of z- variable we have for Berry phase 

WE[n] = -i{2S) / dTzhdrZ (6) 
Jo 

which is apparently independent from the auxiliary co- 
ordinate u. This "simplification" is, however, deceptive 
because the equation n — z^az does not define z un- 
ambiguously. One can always change z(r) — > e*"^^-'z(r) 
with any function Q;(r) satisfying e'"'''^ = e*"'^°). The 
latter is necessary to conserve periodic boundary con- 
ditions = z(0). This "gauge" transformation ob- 
viously does not change the vector n. The change of 
WZNW action (^) under this transformation is 5We — 
i{2S)[a{l3)-a{Q)] = 2(25)277^ where A: is any integer. We 
arrive to the same result: the action We itself is multi- 
valued while e~^^ is single-valued under the condition 
that 2S is an integer numberl 

Let us now briefly discuss topological aspects of Berry 
phase (WZNW term) for a single spin (^. The classi- 
cal unit vector n from (^ takes values on a unit two- 
dimensional sphere - the target space of the problem. 
There exists nontrivial expression (the integrand of (Q)) 
depending on vector n and its derivatives which is (i) 
local (ii) being integrated over two-dimensional sphere 
gives integer values for all possible smooth configurations 
of n defined on S^. Mathematically we say that there ex- 
ists a closed but non-exact 2-form on S"^ (target space) 



III. DOUBLE EXCHANGE MODEL IN THE 
LIMIT OF INFINITE Jh 

Let us now consider the double exchange model (||) 
and take a limit of infinite Hund's constant Jn/t oo. 
We will do it in the Lagrangian formulation. First we use 
the WZNW action ^ to write the Lagrangian (in this 
section we use real time representation) of the double 
exchange model corresponding to the Hamiltonian (^ 



E 



2SzUdtZ, 



- H. 



(7) 



Here /x is a chemical potential which is chosen to guar- 
antee the correct filling factor of electrons and Zi = 
{zii,Zi2Y is a complex vector corresponding to the core 
spin Si at the site i of the lattice. We will make the 
following change of variables in the Lagrangian Ci — 
TpiZi + TpiZi. Here ipi^ipi are spinless fermions, Zi is a 
core spin variable, and Zi is a complex vector orthogonal 
to Zi, namely, Zi = {z*2,—z*if. Physically, -01 is a com- 
ponent of an electron on the site i with a spin parallel to 
the core spin at the same site and ■0^ is the one with a 
spin antiparallel to the core spin. Taking limit Jh — + oo 
corresponds to neglecting "0 component of an electron. 
We substitute Ci = ipi^i info (|^) and omitting constant 
terms obtain 



J2{^S + 4i;i)zjid, 



'tZi 



(8) 



The first term of (||) is a kinetic term for projected 
electrons. Let us take a close look at the second and the 
third terms. We assume that the motion of spins is much 
slower then the motion of electrons. This is justified in 
the limit of large core spins 3> 1 . Then one can think 
of electrons as moving in the background of almost static 
z's. The second term in (||) describes the hopping of these 
electrons with the effective hopping amplitude tl^^^ — 

tzjzj. This amplitude is complex number and has both 
the magnitude and the phase. 

The magnitude of the effective hopping amplitude is 

e//|2 — 

\j is the angle between two core 
spins. In this derivation we used a relation Hi 
and the known identity on Pauli matrices (fafja-ys 



given by ji^/'p ^ = t2(z|z,)(z]z,) = t^^ 
A = t cos ^ where 
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'^SasSpj — Sa/iSjs- The absolute value of the effective 
hopping is proportional to the cosine of the half of the 
angle between core spins It is maximal and equal t 
for parallel core spins and vanishes for antiparallel spins. 

The phase of the effective amplitude t'^j'^ is not de- 
fined unambiguously because of the gauge freedom (see 
the previous section). Consider, however, an electron 
moving along some closed trajectory in a background of 
a smooth configuration of core spins. Then along the 
path we can write z^+i = Zi + 6zi, where 6zi is small, and 
zjzi+i — l + zjSzi. The product of all the effective ampli- 
tudes along the path is Hi^i.'i+i ^ exp(X]i In zj-^i+i) ~ 



; sin Oil 2), where we 



use parameterization = e*'^'(cos6'i/2, e"''* sin0i/2) of 
Zi in terms of Euler angles Oi,(j)i, and ipi. Along the 
closed path J^i ^i'i — 0; so only the second term con- 
tributes to the phase picked up by an electron. It is easy 
to see that this second term is equal to the half of the 
solid angle enclosed by vectors n along the trajectory of 
an electron. We see that (projected) electrons moving 
along some closed trajectory acquire the phase equal to 
the half of solid angle formed by spins along this trajec- 
tory. One can say, therefore, that they "feel" the solid 
angle formed by core spins as a flux of magnetic field. 
This flux is (modulo 27r) a gauge invariant and well de- 
fined quantity. 

Finally, the last term in (||) is a "modified" kinetic term 
for core spins. One might expect that after the averaging 
over electrons this term will describe the dynamics of ef- 
fective spins of magnitude S-V {ip''ip) /2. However, having 
non-integer coefficient in front of Berry phase results in 
an inconsistent theory. 

To make this seeming paradox more precise we average 
over electrons accurately by integrating fermions |l0| out 
of (H). We obtain an effective action for core spins in 
terms of z-variables 



W 



~i Indet |^(5ij(i9t — fJ- + zjidtZi) + tzjz-j 
J dtJ2^Szlidtz,. 



(9) 



Here the expression in the determinant is both an op- 
erator in time and N x N matrix in lattice sites. The 
matrix element zjzj is present only if i,j are nearest 
neighbors. In the adiabatic limit there are two main ef- 
fects produced by an integration over fermionic degrees 
of freedom. First, the energy of the filled (to the filling 
factor x) Fermi sea is added to the action. This leads to 
an effective ferromagnetic interaction between core spins. 
Second, the correction to the Berry phase of core spins 
must reflect the fact that spins of conduction electrons 
are effectively added to core spins. The latter correction 
is the matter of consideration of this paper. 

The variation of (^) over zjidtZi immediately gives 
(25' + i^li/ji)), where {ipjip,) are some fractional num- 
bers which in general depend on the configuration of 
core spins (z^'s). This implies the presence of the term 



{2S + {'ipl^i))zlidtZi in the action (||) and confirms the 
result we expected from the heuristic "averaging" over 
electrons. However, the presence of this term alone ren- 
ders theory inconsistent. 

In the next section we calculate the fermionic deter- 
minant of the type (^ albeit for a simplified double ex- 
change model on two lattice sites. We show that there 
is an additional term present in the effective action for 
core spins. This additional term will make the effective 
theory consistent as it is shown in section 0. 



IV. DOUBLE EXCHANGE MODEL FOR TWO 
CORE SPINS 

The paradox of having fractional spin in the effective 
action for a double exchange model is present already 
for a much simpler double exchange model on two lattice 
sites with a single conduction electron. In addition the 
electron spectrum in such model is discrete which makes 
integration over fermions in adiabatic approximation to 
be a well-defined procedure. Therefore, in this section we 
consider the double exchange model on two lattice sites 
with one conduction electron and show how to resolve 
the issue of fractional spin on this simple example. 

The projected {Jh oo) double exchange model (||) 
on two lattice sites (in imaginary time) 



[>Pl{dr + + (25 + ijlij,)zldrZ, 

1=1,2 

t[i^l{zlz2)i^, 



h.c. 



(10) 



where We = Jo ^'''^e is a Euclidean action. Chemi- 
cal potential /i must be chosen to guarantee that there 
is exactly one conduction electron in the system, i.e., 
J2i=i 2'^i''^-i- ~ 1- ^^'^ '^^^ repeat the heuristic argu- 
ment of the previous section replacing in ths sec- 
ond term of ( p^ by its average value 1/2. One obtains 
two effective spins of magnitude 5-1-1/4 instead of bare 
core spins. Again, we have the problem of having the 
fractional coefficient in the Berry phase term. 
We rewrite the Lagrangian diG) as 



1=1,2 



2Sz]drZi 



+ ¥{^r + l^^ + A)-^. (11) 



Here we defined = ^^"^ 2x2 matrix A 

z\drZi — Ae" 



A 



-Ae"" zldrZ2 



(12) 



where we introduced explicitly the phase and the magni- 
tude of an effective hopping amplitude 



tz\z2 = Ae 



(13) 



The absolute value of an effective hopping amplitude 
A = t cos 6/2, where 6 is an angle between core spins 
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and 5*2. It is equal to zero only for strictly antipar- 
allel core spins 9 = tt. For fixed directions of core spins 
the matrix A becomes constant matrix with eigenvalues 
±A corresponding to the two-level spectrum of an elec- 
tron hopping between two core sites in double exchange 
model with infinite Jh- Two-level spectrum is symmet- 
ric and one can choose the chemical potential /i = so 
that there is only one filled energy level at any given mo- 
ment of time. Then the spinless fermion ijj (electron with 
projected spin) occupies the lowest energy level E = —A 
which is minimal (-Emin = ~t) when core spins are paral- 
lel. Changing the angle between core spins requires the 
energy of the order of t while we are interested in low 
energy dynamics with typical values of energy ~ t/S. 
Therefore, we assume that A does not vanish and is of 
the order of t. 

The Grassman variable c(r) corresponding to the con- 
duction electron must satisfy antiperiodic boundary con- 
ditions in imaginary time r: Ci{f3) — — c(0). We repre- 
sent Ci{T) — 'i(ji{T)zi{T) as a product of spinless fermion 
■01 (t) and spin 1/2 boson Zi{T), which after projec- 
tion Jh — > oo becomes z-boson corresponding to core 
spins. The boundary conditions for 'ipi{(3) = — V'i(O) and 
Zi{(3) = Zi(0) are antiperiodic and periodic respectively. 
Then we have for «'(t) and a(r) (see (|ll|,|l|)) 



*(/3) = -vI'(O), 
a{f3) = a(0) + 2TTk, 



(14) 
(15) 



where k is an integer - the winding number of the phase 



a in imaginary time 27rfc = J^^ dr dr 



The effective action of core spins is given by 

r/3 



eff 



E 

i=l,2 



\2SzldrZ, 



dr 



(16) 



where the effective action of core spins induced by 
fermions 



(17) 



Here the operator dr + A is linear differential operator 
with explicit dependence on time r through the 2x2 
matrix A and we put fi = corresponding to having 
only one conduction electron in the system. To calculate 
the functional trace (|lj) we find the eigenvalues of this 
operator by solving differential equation 



(a, + A)*(T) = a*(t) 



(18) 



with boundary conditions ( |l4| , [l5| ) and calculate the sum 
of logarithms of all eigenvalues A. Matrix (|l^) can be 
rewritten as 



A = -Ao- e-i"" Aa^e^-"" 



(19) 



where 



iAq = z\drZi+ zldrZ2, 



i [z\drZi - zldrZ2 



We would like to solve (|T^) with the matrix A from ( |T9| ) 
in an adiabatic approximation, e.g., assuming that all 
time derivatives are much smaller than the energy scale 
A ^ t. First we "unwind" wavefunction 5" by a unitary 
transformation 



dT(Ao + aa3) X, (20) 



^ = exp 



where d means the derivative of a over r. Then, instead 
of (|l^) the function x satisfies 



(Lo + V)x = Ax, 



(21) 



where we defined "zero order" operator Lq and perturba- 
tion V which is proportional to time derivatives of core 
spins as 

Lo = dr~Aa\ 
V = idacT^, 
1 

Now we can solve ( pl| ) using the perturbation theory in 
V. The boundary conditions of x a-re given by (^4|) and 
(1^) as 

X{(3) = -e*^x(0), 

where Berry phase 

1 ff^ 

where we used that ddr = 27rfc (fc is an integer, see 

(|l|)) and that e*'^'='"' = e"'=. 

We are going to find the eigenvalues of (^l]) using the 
perturbation theory in V. This is justified in an adiabatic 
limit when ao proportional to the rate of change of core 
spins is much smaller than the typical fermionic energy 
A or in other words when the conduction electron is the 
fast degree of freedom compared to core spins. 

The solutions of the zero order equation 

are labelled by an integer number n and by plus or minus 
sign 



\x-:l) 



with eigenvalues 



2/3 



±1 



Ai"2 = ±A+^(2n+l) + l,, 



(22) 



where 
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There are no corrections to these eigenvalues in the first 
order of perturbation theory in V. Therefore, up to the 
second order of perturbation theory we obtain 



A 



— ZTT / 'T 

-Af±A + — 2n+l + - 



(23) 



where we have restored the chemical potential fi which 
should be /i = for the problem with a single conduction 
electron. 

The effective action (RTh is given by 



+ 00 



W^nd^_ J2 [In A, 



In A„- 



(24) 



n— — oo 



To calculate the obviously divergent sum in ( 24 ) we regu- 



larize it by subtracting the effective action for the system 
without conduction electrons. That is we calculate the 
difference between (|2j) calculated at /i = and the one 
calculated at /i — > — oo. 



+ 00 



In- 



A„+(^ = 0) 



In 



^ Xn+ifJ- -oo) 

A„-(Ai = 0) 



-oo 



(25) 



This sum is still logarithmically divergent but one can 
show that the divergent part is a constant not depending 
on physical parameters A and 7. Calculating the sum in 
( p5| ) we obtain 



J 



coshf^+i^'jcoshf^-z^) 

W^nd = - Hm In ( ^ ^. (26) 

cosh(^(^ + *i)cosh(H^-,2) ^ 



In the limit /i — )■ — 00 we obtain 



For the full effective action (Hq) we obtain 



Wind = 17 - In cosh (-^ i- 



/3A 
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In cosh ( ■ h i — 

2 2 



(27) 



This is the final result for the effective action calculated 
at finite "temperature" (3. Notice that the only imag- 
inary part of a Euclidean effective action is the Berry 
phase 7. We have obtained this phase in the first order 
of perturbation theory. However, the same expression is 
also valid in all orders of perturbation theory. It is easy 
to show for our simple model because the higher orders 
of perturbation theory can only give corrections to the 
real part of ( p3|) which does not contribute to the imagi- 
nary part of an Euclidean action. We consider the higher 
orders of perturbation theory using the regular gradient 
expansion in Appendix^. In fact, the imaginary part of 
a Euclidean action has a topological nature and, there- 
fore, can not be changed by small perturbations (see the 
discussion in section ^ . 

In the "zero temperature" limit (3 A ^ 1 we have 



or explicitly 



{Ao + d)dT 



A dr. 



(28) 



The latter result can of course be obtained just by re- 
placing the discrete sum over n in (E5h by an integral. 



dT 



J2 (2S+l)zld,z,+ '-d 



dyn 



AdT 



(29) 
,(30) 

(31) 



The first term of the dynamic part of an effective action 
M^^^" is just the time integral of the energy of the ground 
state of the conduction electron in a fixed instantaneous 
configuration of core spins. The dots denote the higher 
order corrections to (^T]). These corrections are small if 
the background core spins are moving adiabatically and 
they are diverging if derivatives ni , ri2 are bigger than 
A. For such processes an adiabatic approximation fails. 
In particular, if core spins have opposite directions ni — 
—n2 we have A = and an adiabatic approximation is 
not applicable at any rate of the change of core spins. 



BERRY PHASE FOR FRACTIONAL SPINS 



A. Consistency 



The Berry phase term (30) is the main result of this 



paper. We see that the the correction to the Berry phase 
of bare core spins induced by the conduction electron is 
twofold. First, as we expected averaging over fermions in 
^ there is a change of the coefficient in Berry phases 
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of core spins from 25* to 25*+ 1/2. In addition, there is a 
correction ^ J^^^ dr a — ink, where the integer number k is 
a winding number ([l5| ) of the phase a in imaginary time. 
The latter term is a full time derivative and, therefore, 
does not change the classical dynamics of spins. Classi- 
cally, the effective "core" spins behave as if they had the 
value of 5+ 1/4! 

Let us show now that the full Berry phase (^) is 
consistent with quantum mechanics of spins and, there- 
fore, answers the question of how to write down the 
Berry phase for two fractional spins S + 1/4. As we 
have seen at the end of the section || the inconsis- 
tency occurs if the partition function is not invariant un- 
der gauge transformations of z variables. We perform 
the gauge transformation zi^2 — *■ e"i'^('^^zi^2 indepen- 
dently for first and second core spins. In order not to 
change the periodic boundary conditions of 2:1^2 we re- 
quire Aai_2 = cei,2iP) — 011,2(0) = 27r/i 2, where Zi_2 are 
arbitrary integer numbers. Under this transformation the 
three terms of (|30| ) depending on zi, Z2, and a change 
by 27rih{2S + 1/2), 2iTil2{2S + 1/2), and m{li - h) re- 
spectively. It is easy to see that neither of these three 
terms gives an invariant contribution to a partition func- 
tion. E.g., for li = 1 the change of the first term is 
27ri(2S'-|-l/2) which gives —1 being exponentiated. How- 
ever, the total kinetic term ( pO| ) is equal to the sum of 
all these terms and changes under the gauge transforma- 
tion by 2m[{2S + l)/i -I- 2Sl2\ which is a multiple of 2tt 
and, therefore, does not change the weig ht 6"^^= . This 
means that neither of three terms considered makes any 
sense separately from the others [Q. Only their gauge 
invariant (modulo 27r) sum has a physical meaning of the 
combined Berry phase of the system of two core spins 
plus one conduction electron. 

B. Topology and adiabatic approximation 

We have seen in section^ that the existence of WZNW 
term for a single unit vector (single spin) is due to the fact 
that the target space S*^ has a nontrivial second cohomol- 
ogy group H2{S'^) = Z. The addition of WZNW term to 
the action with the coefficient 25 is consistent only if this 
coefficient is quantized (is an integer number). For two 
independent spins the target space becomes S'^ x and 
one obviously can write two independent WZNW terms 
corresponding to H2{S'^ x S'^) = Z x Z . This is appar- 
ently true for two independent spins. However, the result 
is different for the conduction electron induced dynam- 
ics of core spins in the double exchange model. Let us 
discuss why this is so. To obtain the induced dynam- 
ics of core spins we have excluded the charge degree of 
freedom of conduction electron integrating it out in an 
adiabatic approximation. This approximation, however, 
breaks when the the rate of the change of core spins be- 
comes comparable to the energy scale 2A given by the 
levels of conduction electron in the background of time 
independent core spins. There is a special point though 




FIG. 3: The contraction x S^\S^ is shown as a one- 

dimensional illustration of G 5*^ . a) The torus space x 
is represented as a square with the opposite sides identified 
according to the direction of arrows. The diagonal shows the 
points given by ni = — n2. b) The space x S^\S^ — > 
is obtained by the removal of the diagonal, c) One pair of 
opposite sides is glued, d) The second pair is glued making a 
cylinder, e) The cylinder is contracted onto the unit circle. 

when core spins are opposite to each other. At this point 
A = and an adiabatic approximation breaks no matter 
how slow the motion of core spins is. 

This is quite obvious from the physical point of view. 
Consider two time independent core spins opposite to 
each other. The effective hopping amplitude of a con- 
duction electron is ~ cos 61/2 — and one has two dif- 
ferent configurations. The electron can be either on the 
left lattice site or on the right one (see Figure ||) and 
the matrix element of the hopping term is zero between 
those states. This means that the ground state of the 
system is doubly degenerate. The charge degree of free- 
dom (left or right position of the electron) is essential for 
this particular configuration of core spins and can not be 
integrated out. In our derivation we used an adiabatic 
approximation and neglected such configurations. This 
is justified in the limit 5^1 for ferromagnetically in- 
teracting core spins. However, this approximation dras- 
tically changes the topology of the system. Essentially, 
we prohibited the configurations ni = —n2. The relevant 
target space changes from x for two independent 
spins to G = 5^ X S^\S^. The latter expression de- 
notes the direct product S'^ x with points ni = — n2 
(topologically equivalent to one more S'^) removed. The 
resulting target space G is topologically equivalent to 
S^. To prove it one can show that it can be con- 
tracted along itself onto 5^. We illustrate the similar 
phenomenon on the corresponding one-dimensional ex- 
ample in Figure ^. Therefore, the relevant target space 
G ^ S^. To check for allowed WZNW terms we calculate 
H2{G) = H2{S^) = Z. There is only one WZNW term 
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allowed for the chosen target space! 

The existence of only one allowed topological term ex- 
plains why only the full Berry phase (|o|) is well defined. 
The most general form (in CP^ representation) of the 
topological term for two spins is given by 



Dmitri Ivanov is deeply appreciated. A. G. A. is Alfred 
P. Sloan Research Fellow. The work of Ar. A. is sup- 
ported by DR Project 200153 at Los Alamos National 
Laboratory. 



dr 



pizldrZi + p2zldrZ2 + i£.d ■ (32) 



Coefficients pi,2,S, satisfy quantization conditions (or 
gauge invariance conditions). Namely, pi — ^ and P2 + C 
must be integers. We see that for the double exchange 
model of two core spins and a single conduction electron 
we have (|o|) which corresponds to particular integer val- 
ues pi — £, = 2S and p2 + £, = 2S + 1. If by some reason 
the conduction electron spends unequal time at two sites 
(e.g., the local magnetic field acting on electron spin on 
one of sites is added) we have an asymmetric situation 
Pi P2- However, the combination pi— ^ and P2+C must 
stay integer to guarantee the consistency of the theory. 



APPENDIX A: PATH INTEGRAL FOR SINGLE 
SPIN AND FERMIONIC DETERMINANT 



In this section for the sake of completeness we rederive 
the well-known path integral representation for the quan- 
tum mechanics of a single spin. We reduce the problem 
to a calculation of a simple fermionic determinant. The 
method we use differs from the commonly used coherent 
state representation (see e.g., 0) for the path integral. It 
is conceptually simpler but involves the gradient expan- 
sion calculation of a fermionic determinant. It turns out 
that the similar determinant is needed for our analysis of 
a double exchange model for two core spins. 



VI. CONCLUSION 

In this paper we considered the infinite Hund's cou- 
pling limit of the double exchange model with a given 
concentration x of conduction electrons. In the limit of 
large core spins S" ^ 1 one can think of electrons and 
core spins as of fast and slow degrees of freedom respec- 
tively. The superficial averaging over the fast motion of 
electrons leads to an inconsistent theory with effective 
fractional spins. Using the example of the double ex- 
change model on two lattice sites we showed that in ad- 
dition to the change of the value of the effective core spin 
on site from S to S + x/2 'A new topological term arises 
in the effective action. With this term the exponent of 
the effective action becomes single-valued and well de- 
fined. We showed that the full Berry phase (|30| ) reflects 
the change of topology of the system which occurs due to 
the adiabaticity condition — the condition justifying the 
averaging over the charge degrees of freedom. The gen- 
eralization of our results to the double exchange model 
for many spins is relatively straightforward. The only 
major difference with the case of two spins is that the 
average occupation number of a given lattice site by con- 
duction electrons depends on the configuration of core 
spins. This makes the magnitude of a "fractional" spin 
on the given site a fluctuating quantity with correspon- 
dent modifications of Berry phase terms. The latter ef- 
fect is not expected to be strong in the ferromagnetically 
ordered state. 
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1. Path integral representation for spin S 

It is easy to check that the representation S = 
5V'i'?a/3^/3 with constraint ipl^4'a = 1 is a faithful rep- 
resentation of (||) with S = 1/2. Here a = (cr^, cr^, cr^) 
is a set of Pauli matrices and ipajfp}) with a, /3 = 1,2 
are annihilation (creation) Fermi operators respectively 
which satisfy the following (anti) commutation relations 

{?Aq,i/'/3} = = and {ipa,^l} = ^ap- In fact, 

this representation reflects how spin operators appear in 
Nature with being an annihilation operator of the phys- 
ical electron. 

We consider the following Hamiltonian H = —2AnS, 
where S — ^ip^aip is a spin of a fermion and A is a 
strength of the coupling of this spin with the unit vector 
n. The correspondent action is 

S = J dt^^ [idt + An] V^, (Al) 

where we introduced a notation n = na. One can think of 
2An as of the magnetic field acting on the fermionic spin 
S. In the limit of an infinite coupling A ^ oo (infinite 
magnetic field) the fermionic spin S will be "frozen" along 
the classical vector n so that n in the limit of infinite A 
becomes a direction of quantization axis of a quantum 
spin of a fermion. If now we change the direction of the 
classical vector n with time, the fermionic spin will follow 
the direction of n. Therefore, we expect that in the limit 
of an infinite A an effective dynamics of classical vector 
n induced by a fermion will be that of the quantum spin 
1/2. 

The path integral representation for spin 1/2 is then 
given by 

Z ^ [ £lne'^["I, (A2) 
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where W[n] = liniA^cx) W^^f[n] is a classical action cor- 
responding to a quantum spin 1/2. This classical action 
is a functional of the trajectory of the time dependent 
unit vector n{t). The effective action VFe//[ri] is defined 
as 



Dil}D'4>'^ e 



t i J dttp''[idt+An]tp 



(A3) 



The functional inte gra tion in (A3) is performed over 
Grassman variables |10| ip,^^ - classical analogs of 
fermionic annihilation (creation) operators. Integrating 
over fermions we obtain 



W 



eff! 



- i In Det [idt + An] , 



(A4) 



To generalize the result to an arbitrary spin S we can 
consider the same fermionic theory ( |A3| ) but with 25 in- 
dependent species of fermions coupled to the same vector 
n. In this case all s pins 1/2 of fermions will be aligned 
along n so that (A2) will correspond to quantum spin S. 
Integration over all fermionic species is independent and 
we have for an effective action 



Wf^irt] = -i2S In Bet [idt + At 



(A5) 



The only problem which is left is technical. It is to cal- 
culate the fermionic determinant (A5) in the limit of an 
infinite A. 



2. Shortcut 

One can obtain the result without (almost) any calcu- 
lations. Let us introduce the complex, two-component, 
normalized vector z* = (21,22), z'^z = 1 so that it is an 
eigenvector of h with the eigenvalue 1: nz = z. For a 
time independent fi this vector is a ground state of our 
Hamiltonian Hz = —Ahz = —Az with the ground state 
energy —A. The other eigenstate has an energy -|-A and 
in the limit of an infinite A is never excited by a motion 
of the unit vector ft. Therefore, in the limit of an infi- 
nite A we can use ipa — ZaX with x ^ spinlcss (projected) 
fermion. Substituting this into the action ( |Al[ ) we obtain 



S 



dt [xHidt + A)x + xh{z^idtz)\ 



The variable x does not have any dynamics in the limit 
of an infinite A. We average over x using x^X = 1 ^^'^ 
obtain W = J dtzHdtz for spin 1/2. In case of spin 
S one has 25* species of fermions, x^X — 25* and W = 
2S J dt zHdtz. This is the desired result for the action for 
a single spin S. However, the main point of this paper is 
to show that the "naive" averaging over fermions instead 
of the careful determinant calculation can result in error. 
In this case one can show that due to topological reasons 
the calculation we performed is a legitimate one. To add 
even more weight to this statement we procee d w ith a 
formal calculation of the fermionic determinant ( A5 ) and 
show that it gives the same result for a single spin action. 



3. Determinant calculation 

Bef ore going to the calculation of fermionic determi- 
nant (A5) we will perform Wick rotation of time and 
go to an imaginary time representation. Imaginary time 
representation is slightly more convenient for calculations 
for reasons that will be clear later. The change to an 
imaginary time changes quantum amplitude into a 
"Boltzmann" factor e^^^, where subscript "E" stands 
for "Euclidean" . In imaginary time t — it we have in- 



stead of ( |A5D 



W^^^[n] = -25" In Det [a^ - An] 
with the partition function for a single spin 



Z 



Dn e-^^["l 



(A6) 



(A7) 



where W£;[n] ~ \imA-*ooW^ [n]. We impose periodic 
boundary conditions (in imaginary time) on n(r) so that 
n(/3) — n(0) and /3 is a maximal span of an imaginary 
time. With these boundary conditions we can think of 
time as of the one-dimensional circle with identification of 
points T = /? with r = 0. Physical results will be obtained 
in the limit (3 —^ 00 after analytical continuation back to 
real time t = —ir. 

Using the method presented in Appendix |b| one can 
calculate the effective action ( |A6| ) in adiabatic approx- 
imation with an arbitrary accuracy. However, we are 
interested in the limit A ^ 00. In this limit all terms 
of gradient expansion in 1/A except for the first (topo- 
logical) one will vanish. Here we show how to obtain 
the most important (to pological) A-independent part of 
the effective action (A6) which is the exact result for the 
action for a single spin. 

Determinant of an operator does not change if one ap- 
plies a unitary transformation to it. This is, however, 
a potentially dangerous procedure in the context of gra- 
dient expansion as one has to make sure the boundary 
conditions for electrons do not change in the process of 
such transformation. It turns out that this is not an issue 
for the single spin problem. However, e.g., for the dou- 
ble exchange problem the boundary condition of fermions 
do change in a non-trivial way under the unitary rotation 
(see section IV). 

Let us use the transformation U such that U^^fiU = 
a^. This transformation U is easily found using the so 
called CP^ representation. In this representation we in- 
troduce a two-component complex vector z with the con- 
straint z^ z — 1 such that n = z^'az. Then the matrix U 
can be written as 



U 



_ I Zl Z2 

zi -zl 



We write the equation (A6) in the form 



2S'Trln dr + A-Aa^ 



(A8) 



(A9) 
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where A = U~^drU and Tr means both functional trace 
and the trace over a-matrices. 

The matrix elements of A are proportional to the first 
time derivative of spin variables and are small incom- 
parison to A. We can expand the logarithm in (A9) in 
powers of A. We keep only the first A dependent term 
as the other terms will have additional powers of A in 
denominators, (one can obtain those terms in a regular 
way explained in Appendix ^) 



Wl."[n] = 25 Tr 



dr + Acr3 



A 



duj 
2^ 



■tr- 



A2 



A{uj = 0) 



24 



/3 



dTtva'^A 



/3 



i2S I dr zHdrZ. 



(AlO) 



It is instructive to rewrite ( AlO ) in terms of the 
original variables n. We parameterize complex vector 
z* = e''^(cos |, e"^ sin |). It is easy to see that the an- 
gles 9 and (f) are the polar and azimuthal angles cor- 
responding to the direction of the unit vector n. The 
phase tp can be chosen arbitrarily as it does not af- 
fect n. Up to the full time derivative of tp we obtain 
-J^dTzUdrZ = 5/(f dT(l - cos 61)0 = n/2, where n 
is a solid angle enclosed by the unit vector n during its 
time evolution. One can check that this solid angle can 
be rewritten in the form (^ which is explicitly rotation- 
ally invariant. 



where 



SA = iSAa/2 + SAa. 



Let us now calculate the functional trace in the basis of 
plane waves 



tr 



2tt J dr + A 



(Bl) 



Here the integral over time t is the calculation of the 
diagonal matrix element {uj\5DD~^\uj) and integration 
over oj is the calculation of the functional trace. The 
residual trace tr is the trace in the space of 2 x 2 matrices. 
"Pulling" e*'^'^ from right to the left we have 



tr 



^ fdrSA ^ (B2) 

2n J iuj + dr + A 



In the latter expression the partial derivative dr is as- 
sumed to act "to the right" . This means that when we 
expand 



tr / dt 



duj 



SAG [1 - drG + drGdrG 



= j dt[SLo 



2tt 

6L1 + 6L2 



(B3) 



all derivatives are assumed to act on all terms to the 
right. E.g., drGdrG = GG + GG. We introduced 



G 



1 



-i{uj + Ao/2) + Act 
^ " (co + ^o/2)2 + |A|2 



(B4) 



APPENDIX B: CALCULATION OF FERMIONIC 
DETERMINANT FOR DOUBLE EXCHANGE 
MODEL 

In this appendix we calculate the effective action ( p^ ) 
which can be written as 



W, 



ind 



Tr InD, 



and 



SLi 
SL2 



tr 



tr 



tr 



diu 
2^ 



SAG, 



^SAGdrG, 
2tt 



duj 
2^ 



SAGdrGdrG. 



(B5) 



where the differential operator D is defined as 
D^dr+A 



with matrix A given by (g^|jl9|). We rewrite 

A = iAf)/2 + iao(T^ + Aict^ + A2(t^ = iAn/2 + Act, 

where Ai — — A cos a, A2 = — A sin a, and A3 — mo. 

We use an adiabatic approximation using the (gradi- 
ent) expansion in 1/A. The operator D contains the 
magnitude of the effective hopping A which is slowly 
varying but it is not small itself. To avoid this difficulty 
we calculate instead of effective action its variation and 
then restore the result from this variation 



SWi^"^ = - TtSDD-^ = - Tr 



SA- 



1 



dr + A 



We notice that lAP = A^ 



A^ because ao ^ A by 



adiabatic approximation (it contains time derivative of 
core spins). Therefore, the denominator of (B4) is never 
zero and integrals (B5) are not singular. The expansion 
(B3) is essentially the gradient expansion, i.e., the expan- 
sion in the number of time derivatives. One can see that 
Li contains at least one time derivative, L2 - at least two 
etc. 

The calculation of SLo is straightforward and gives 



or removing variation 



iAo/2-\A\ 



(B6) 



2 2A ^ ' 
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To obtain the latter expression we expanded up to the 
second order in time derivatives using | A| = ^ — Oq k, 

Let us now prove that the whole dependence of the 
full result ( p3| ) on comes from Lq. Varying the action 
(|ll|) over Aq we obtain |(vl/tvli) = i. This statement is 
exact because the total number of electrons in the system 
is fixed (^^vj;^ = \ por the effective action we have, 
sw^^" ^ However, from (p^"^ — 



therefore, 



we have 44^ = 



and all higher order terms 



h which gives 

Li, L2 etc. do not depend on Aq. This simplifies further 



calculations allowing to use 



G = 



lUJ - 



A 



Act 



Ao=0 



|A|2 



for all higher order terms. 

As a result of calculations one obtains from 



Li = ia/2 which gives half of the solid angle after time 
integrati on [p^ . Calculating corrections due to nonzero 
flQ from (B8) we obtain 



it. la^a 
2 I lAI 



2""" 2A' 



where we expanded the latter expression up to the terms 
of the second order in time derivatives. 



The calculation of L2 in (B5) is straightforward and 
gives 



8IAI I "lAI / ~ 8A- 



(BIO) 



Finally, combining terms (B^,B£. B10| ) we obtain an 
effective action up to the second order terms 




(B8) 



This is identical in form to a Berry phase of a single spin 
1/2 given by (||) with identification n — A/|A|. There- 
fore, Li is given by formula analogous to the Berry phase 
for a single spin and is to be interpreted as a half of 
solid angle swept by vector A during its time evolution. 
Assume for a moment now that ao — 0. Then the vector 
A/|A| = (— cosa, — sina, 0) moves along the equator of 
unit sphere and covers the solid angle — a{f3) — a(0) 
during its motion. In this approximation we obtain 



{Ao + d) 



(Bll) 



It is intere sting to not ice that although neither of con- 
tributions (B7 Bgj . BlC ) is gauge invariant, the full action 



(Bll ) is g auge invariant. This means that the final re- 



sult (Bll) can be rewritten in terms of two unit vectors 
ni and n2 which represent physical degrees of freedom - 
core spins of double exchange model (M. 

The method of gradient expansion used in this ap- 
pendix can of course also give higher order corrections 

to (pnl). 
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